opt=odeset("OutputFcn","odeplot");
ode45(@rhs,[0,50],-10,opt);
print(['doc/dotx=-satx','.eps'],'-deps','-r600');% save eps with 600dpi
print(['doc/dotx=-satx_epsc','.eps'],'-depsc','-r600');% save eps with 600dpi
print(['doc/dotx=-satx_eps2','.eps'],'-deps2','-r600');% save eps with 600dpi
print(['doc/dotx=-satx_epsc2','.eps'],'-depsc2','-r600');% save eps with 600dpi


function dxdt=rhs(t,x)
    dxdt=-sat(x,-1,1)+1.1;
end